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A nonstandard method for calculating the nuclear 2vflf3 - decay amplitude is proposed. The 
method is based on the explicit use of those approximate symmetries of a nuclear hamiltonian, 
which correspond to the operators of allowed (3 - transitions. Within the framework of the proposed 
- - . method the mentioned amplitude is calculated for a wide range of nuclei. The model parameters 
Q\ • used in calculations are taken from independent data. Calculated 2v(3(3 half-lifes are compared 



with known experimental data. 
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I. An analysis of a great body of known experimental data (see e.g. ref. |T| ) allows one to conclude 
that the 2v(3(3 - decay rate is hindered as compared with the rates evaluated within the independent- 
^ ■ quasiparticle approximation. The hindrance is caused by existence of the isobaric analogue and 
^ ; Gamow- Teller nuclear collective states (IAS and GTS) with relatively large excitation energy. These 
— i states each exhaust the main part of the Fermi and Gamow- Teller strength, respectively. To describe 

in ■ 

O ' the hindrance the different versions of the quasiparticle random phase approximation (QRPA) are 
00 ' 

Q | used for calculations of the nuclear 2vf}(3 - decay amplitude. It has been found (see e.g. refs. |2|, ^]) 
! that the calculated amplitudes are very sensitive to the particle-particle interaction strength. For 
^ this reason the predictable power of the QRPA current versions is poor. As reported recently, the 
\ similar instability of the calculated Oz//3/3 - decay amplitude has been also found || |J. We believe 
^5 ! that the mentioned difficulty caused by inconsistency of the QRPA current versions can be presently 
^ overcome by the explicit use of the approximate isospin and spin-isospin symmetries of a nuclear 
\ hamiltonian. Early attempts along this line have been undertaken in refs. [0, [|. 

It is well-known that the Fermi-type excitations can be correctly described only with consideration 
for the isospin symmetry of a nuclear hamiltonian. Such a consideration allows one to reproduce 
in model calculations the experimental fact, that IAS exhausts more than 95% the Fermi sum rule 
N — Z, a.s well as to evaluate correctly the small Fermi-strength of the low-energy IAS satellites 
populated virtually in the 2vf3f3 - decay process. If the isospin symmetry were exact, the IAS would 
exhaust 100% of the Fermi strength and the Fermi nuclear amplitude for the 2vf3(3 - decay to the 
ground or low-excited states of a product nucleus would be equal to zero. The above statements 
allow one to describe the Fermi-type excitations within the framework of the perturbation theory 
with respect to the variable part of the nuclear mean Coulomb field. This latter is the main source 
of the isospin-symmetry violation in intermediate and heavy mass nuclei (see e.g. ref. Q and refs. 
therein) . 
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The similar situation takes place at the description Gamow- Teller (GT) excitations. It is experi- 
mentally established that in intermediate and heavy mass nuclei the GTS exhausts about 70% - 80% 
sum rule 3e^(N — Z), where e q ~ 0.8 is "a effective charge" describing, in particular, the renormaliza- 
tion of the axial constant of the weak interaction in nuclei. The mentioned experimental fact can be 
considered as an evidence of approximate conservation of the appropriate spin-isospin symmetry of 



a nuclear hamiltonian (see in this connection ref. fllP]). This symmetry should be taken into consid- 



eration for correct description of both the GTS low-energy satellites and the nuclear 2vf3(3 - decay 
amplitude. If the mentioned symmetry were exact, the amplitude would be equal to zero. Because 
the 2vf3f3 -decay hindrance is connected with formation of the collective GTS, the hindrance can 
be also considered as a result of partial conservation of the appropriate spin-isospin symmetry. 

In the presented work a method for calculating the GT nuclear 2vf3f3 - decay amplitude is pro- 
posed. The method is based on the explicit consideration for approximate conservation of the sym- 
metries corresponding to operators of the allowed Fermi and GT j3 - transitions. The 2vf3~ (3~ and 
2vf3 + f3 + (including the electron capture) half-lifes for a wide range of nuclei are calculated within the 
framework of the simplest version of the proposed method. The model parameters used in calculation 
are taken from independent data. Calculation results are compared with known experimental data. 

2. To clarify the main point of the proposed method we derive the basic formula for nuclear 
2v(5(5 - decay amplitude by the same way for both the Fermi and GT amplitudes. At first we 
consider the 2v(3~/3~ - decay. Let \i) and |/) be the wave functions of the ground state, respectively, 
of a double-even parent-nucleus (N + 2, Z) with energy Ei and of a product-nucleus (N, Z + 2) 
with the energy Ef, be the operators of allowed Fermi (G^ = J2 a T i a )) or Gamow- Teller 
(G(±) = E a v{a)T {±) {a)) (3 - transitions. 

We start with the zero approximation, in which symmetries corresponding to operators G^ are 
exact. In this approximation the states \G) (IAS and GTS) of an isobaric nucleus (N + 1, Z + 1) : 

(i) each exhaust 100% of the corresponding sum rule, so that 

\G) = k G G«|z>, G«|i) = G«|/> = 0; k G = (z|[G« G(-)]l^)~ 1/2 ! (1) 

(ii) have the same energy Eg] (hi) are orthogonal to other states \A) of the same isobaric nucleus: 

(G\A) = 0- G i+) \A} = 0. (2) 

For the further analysis of the nuclear 2vf3f3 - decay amplitudes it is necessary to consider also the 
following excited states of a product-nucleus: 

\GA) = k A G^\A), E GA = E G -E i + E A , k A = (A\[G< + \ G^]\A)~W , (3) 
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\DG) = k D GHGH|z), E dg = E G -E l + E G , k D = (i\G^GWG^G^\i)- 1 / 2 . (4) 

Low-energy states \A) are called the anti-IAS or anti-GTS, whereas states \DG) are called the double- 
IAS or double-GTS. 

Let U be the nuclear mean field, so that U = £ a U(r a ) and [U, G^) = U G ' ] ^ 0. The part of 
the nuclear mean field U G = £ a U G (r a ), which leads to nonzero operator U G , is the main source of 
violating the symmetry connected with operator G^~\ Other sources are less essential and briefly 
considered below. In the lowest order in the U G the energy of state (1) is equal to 

E Q -Ei = k 2 G (i\[G W , U G - ] ]\i) = A G . (5) 
Let us turn to the expression for the nuclear 2pj3~ (3~ - decay amplitude (see e.g. refs. fll], HU): 



M G = Y.(f\G { - ) \S)(S\G^\i ) u- s \ (6) 
s 

where the sum is taken over intermediate states of an isobaric nucleus (N + 1, Z + 1); oj$ — E$ — 
(Ei + Ef)/2 is the excitation energy of the intermediate state. The 2v(3 + (3 + - decay amplitude is 
equal to conjugate value M G . In the zero approximation the 2vf3f3 - decay is doubly forbidden 
(Mq = 0), as it follows from eqs. (l)-(6) after replacing \S) — > \A), \G). Mixing of nuclear states due 
to the single-particle field U G results in nonzero value of amplitude M G . We keep those corrections 
to wave functions of isobaric nuclei, which provide the contribution to the amplitude M G of the same 
order perturbation with respect to U G : 

\G) -> \G) + Y, a AAA); \A) \A) + a ijA \G) ; a Ajl = -a* A ; 

A 

(7) 

\f)^\f) + j:a A j\GA) + p itf \DG) . 

A 

With the help of eqs. (l)-(4) the first-order mixing amplitudes a and the second-order amplitude f3 
can be expressed in terms of matrix elements of the single-particle charge-exchange field U^: 

OA _ Wm _ k (4ft) n * _ (f\U\GA) _ u iMPlA) 

(8) 

v (f\U\GA)(GA\U\DG) ^ (f\U^\A)(A\U^\i) 
PiJ ~^ (E G a — Ef)(EuQ — Ef) Kd ^ (u g + u a )u g 

After substitution of the eqs. ©,© in eq. (||) with the use of eqs. (1)— (4) we get: 
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M G = ^Y.{f\UP\A){A\U^\^ . 

A 

This equation can be presented in the form, which is similar to the starting expression (6): 

M G = Uo 2 j:(f\vt ) \S)(S\vt ) \i)uJs 1 , (9) 
s 

where V G ~ ] = U [ G ] - A G G^\ A G is determined by eq. (5). Indeed, the equality (G\V G = and 
eqs. (1),(2) allow one to include formally the state \G) in sum (^). 

The fundamental character of formula @ is explained by the following, (i) In this formula 
the approximate symmetries of a nuclear hamiltonian are explicitly taken into account, (ii) The 
hindrance of the strength of intermediate (3 - transitions due to formation of the collective states 
(IAS and GTS) is correctly taken into consideration in rather general form, (iii) Inasmuch as the 
main part of polarization effects caused by a quasiparticle interaction has been taken into account, 
rather simple nuclear models can be used for calculating the intermediate-state wave functions and 
energies in eq. (|9p. 

3. As the first step in evaluation of the 2v(3(3 - decay amplitude according to eq. ([5J) we 
use: (i) realistic phenomenological mean nuclear field with consideration for partial self-consistency 
conditions; (ii) the BCS model for nuclei with strong nucleon pairing; (iii) the pair-vibration model 
for description of the " magic ± two nucleon" subsystem. 

We use the following phenomenological quantities involved in the model hamiltonian. The 
isoscalar part of the nuclear mean field Uq(t) is the sum of central and spin-orbital parts: 

U (r) = U {r) + U so (r); U so {v) = U so (r)(ls) . (10) 
As an isovector part of the particle-hole interaction the Landau-Migdal forces are used: 



^(n,r 2 ) = {F' + G'a l a 2 )TxT 2 5{v l -v 2 ) , (11) 

where F' and G' are phenomenological constants. Other parts of the nuclear mean field are calculated 
by the self-consistent way. As a result of the isospin symmetry of the nuclear hamiltonian, the 
self-consistency condition, which relates the isovector part of nuclear mean field Usy = |r^v(r) 
and neutron-excess density p(r) = p n (r) — p p (r) via the interaction (|i~lf), takes place (see e.g. ref. 



131): v(r) = 2F'p(r). The nuclear mean Coulomb field Uc(r) is calculated within the Hartree 



approximation via the proton density p p . 

As mentioned above, the mean Coulomb field U c = E a U c {r a ){l ~ r {3) (a))/2 is the mam source 
of the isospin-symmetry violation in intermediate and heavy mass nuclei. In particular, this field 
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results in the shift of the IAS energy from energy Ei. According to eq. (5) the shift equals the 
Coulomb displacement energy Ac-: 

E IAS - Ei = A c = (N + 2 - Z)- 1 J U c {r)p{r)dv. (12) 
Here p(r) is the above-mentioned neutron-excess density in the ground state of the parent nucleus: 

ftr) = ^ (e^W( 2 > + *H -E^W( 2 ^ + iKJ , (i3) 

where R\(r) are the radial single- neutron (A = v) and single-proton (A = n) wave functions; A = 
n r ,l,j is the set of the single-particle quantum numbers; ri\ are occupation numbers satisfying the 
equations: 

E(2> + IK = N + 2, J2( 2 ^ + !K = Z. (14) 

V IT 

For nuclei with strong pairing in any nucleon subsystem the occupation numbers are n\ = v\ = 1— u\, 
where v\ and u\ are the Bogoliubov-transformation coefficients. In the case that a nucleon subsystem 
is "magic ± two nucleons" one the occupation factors are: n\ = n™ + c|(l — 2n™), where n™ are the 
occupations numbers for the magic subsystem; c\ are the coefficients determining the pair-vibration 
wave function and satisfying the normalization condition: J2\ c\(l — 2n™)(2j\ + 1) = ±2. 

In view of a high degree of the isospin-symmetry conservation in nuclei we omit the analysis 
of the relatively small Fermi- 2v(3(3 - decay amplitude and turn to the GT amplitude Mqt- The 
mean Coulomb field results in the same (in the lowest order of the perturbation theory) energy 
shift Ac (12) for the both IAS and GTS. However, the main source of violation of the spin-isospin 
symmetry corresponding to operator = J2a ^a r i ^ is the spin-orbital part of nuclear mean field 
(10): Uso — J2 a Uso( r a)- In particular, this field results in the energy shift A so of the GTS from 
the IAS energy. According to eq. (5) we have: 

Aso = ^ N + 2 - Z) ^ G{+) ^ Uso ^ G{ ^ = - 3(N+\-Z) ^ Uso ^ 

(15) 

(i\Uso\i) = E (MUso(r)\X)(ls) jxlx (2j x + l)n, . 

X=TT,V 

Here (n\Uso( r )\ l/ ) is the radial matrix element; n\ are occupation numbers satisfying eqs. (14). 

In the case of GT excitations the mixing (spin-orbital) field in eq. (9) has the form: Vq ^ = 
H a Vso{ a ) T i~\ where Vso(a) = Uso( r a) [Is, er] — A so <j. Calculation of amplitude M GT according to 
eq. @ within the framework of the BCS model results in the expression: 
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M C rr = elLJ c ^ 3 ^^2l v + l){j v -j u )('ir\Uso(r)\v) - A S o(ir\v)J {ix^a^vfu^v^v^^ . (16) 
Here (it\v) = j " R n (r)R v (r)r 2 dr is the overlap integral; (7r||cr||v) is the reduced matrix element; 



u)^ v — £ n + E v is the excitation energy of the two-quasiparticle state, £\ = J (e\ — fi) 2 + A 2 is the 
single-quasiparticle energy for subsystem with strong nucleon pairing, e\ is the energy of single- 
particle level, A is the energy gap. If nucleon subsystem is magic in final (initial) state and "magic ± 
two nucleon" in initial (final) state the following replacements in eq. (16) should be made: u\v \ — > c\ 
and £\ = |A+eA — £i| for particle pair- vibrations or S\ = |A+e — e A | for hole pair-vibrations, where €\ 
(eo) is the energy of the first empty (last filled) single-particle level in the magic subsystem; 2A is the 
pair- vibration state energy, which coincides with the pair energy P for the subsystem "magic ± two 
nucleons" . The hindrance factor h G T, which can be considered as a perturbation theory parameter, 
is estimated as follows: h GT = M GT / M GT , where amplitude M GT is evaluated accordind to eq. (^) 
within the framework of the BCS model without consideration for the spin-isospin symmetry: 

M GT = e l Yl ( n I V ) 2 ( n II a \\ V ) 2 Ut T V u U v V ij U)~1 . (17) 

In the conclusion of this section we briefly consider other sources of the spin-isospin symmetry 
violation of the nuclear hamiltonian. The inequality F' ^ G' in eq. (|ITD results in appearance of 
the mixing field Vt ] = a v(r a )T( ^(a), where 5 = (G' — F')/F', v(r) is the symmetry potential 
considered previously. In view of smooth radial dependence of v(r) (as well as in the case of mean 
Coulomb field Uc) the contribution of field Vg to the GTS energy should be only taken into account 
1(J. By analogy to eq. (12) we have: 



A s = 5(N + 2 - Z)- 1 J v(r)p(r)dr . 

Thus, the GTS excitation energy in eq. ([[£]) equals uj G ts — Ac + A so + As + \{Ei — Ef). 

Another source of the spin-isospin symmetry violation is the particle-particle interaction and, in 
particular, the nucleon pairing. The relative contribution of this source to the violation intensity 
can be estimated as the ratio A njP / A/ s ( A/ s is the mean energy of the spin-orbit splitting near the 
Fermi energy), which is rather small (< 20%) for intermediate and heavy mass nuclei. Thus, we do 
not expect of the essential dependence of amplitude Mqt calculated within the framework of the 
proposed method on the particle-particle interaction strength in contrast to the standard methods 
(see e.g. §). 
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4. Let us turn to the choice of model parameters and to calculation results. The parametriza- 
tion and phenomenological parameters of the isoscalar part of the nuclear mean field are given in 
detail in ref. |Tf|. The intensity F' = 300 MeV fm 3 in eq. (JO]) is chosen to describe the ex- 
perimental neutron and proton bound-energy difference for nuclei with rather large neutron excess 



48 Ca, 68 Ni, 132 Sn, 208 Pb [14|, for which the difference is mainly determined by the symmetry poten- 
tial and mean Coulomb field. The intensity G' = 255 MeV fm 3 in eq. (ITT]) is chosen according 



to ref. fl~7j] . The parameters /i n)P and A„ iP , which are used for calculation of the Bogoliubov- 
transformation coefficients, were found for each subsystem with strong nucleon pairing so that to 
satisfy eqs. (14) and to describe the experimental pairing energies P = 2 minx £^ taken according to 
ref. 0]. The pair- vibration state energies 2 A = P are calculated on the experimental pairing ener- 



gies [14]. Coefficients c\ determining the pair- vibration state wave function for a nucleon subsystem 



"magic + or - two nucleons", are calculated according to equations c\ ~ (1 — 2n\)/(e\ — t\ + 2A) 
or c\ ~ (1 — 2n\)/(e — t\ + 2A) with taking into account the normalization conditions given in 
Sect. 3 ||15|| . For all considered nuclei the single-particle basis including all bound states as well as 
the quasi-bound states up to ~ 5MeV is used. If the nucleon subsystem in a parent nucleus is 
"magic ± two nucleons", and in the product-nucleus is "magic ± four nucleons", the calculations 
were performed within the framework of the BCS model with the use of u\, v\ factors calculated for 
a "magic ± four nucleons" subsystem. 

Calculation results and corresponding experimental data for a wide range of nuclei are given in 
Table 1. The M GT amplitudes calculated by eq. ( |T6| ) are given (in m" 1 , m e is the electron mass) 
in column 5. The M GT amplitudes were calculated by eq. (|I7]) with the use of the same model 
parameters and the same number of basic states. The calculated hindrance factors her are given in 
column 6. Half-lifes Tyif calculated by formula (T^) -1 = G2v\Mgt\ 2 are given in column 8. The 
lepton factors G2 V taken from refs. ]l|, [TTJ, [12| are also given (column 3). We calculated also the 
2vf3 + f3 + half-lifes (including the electron capture) for those nuclei, for which appearance of relevant 
experimental data is expected. 

Two main conclusions follow from the data given in Table 1. (i) The calculated hindrance factors 
h G T are found to be sufficiently small to justify the use of the perturbation theory with respect to 
the spin-orbital part of nuclear mean field, (ii) Except for the Te isotopes, the calculated amplitudes 
Mqt within the factor 2-3 are in agreement with corresponding amplitudes M^j? deduced from 
experimental data. 

The following step in the analysis of the 2vf3f3 - decay amplitudes by the proposed method 
can be the evaluation of amplitude (9) with the use of the intermediate-state energies and wave 
functions calculated within the QRPA. As mentioned above, the main reason for the hindrance of 
the 2u(3j3 - decay amplitude, which is caused by existence of the GTS, has been taken into account 
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in eq. (9). That is why we do not expect of essential change of amplitudes Mqt as compared with 
calculated ones. A consideration for the particle- hole interaction (for instance, in form (pTT|) ) will 
cause to some redistribution of the small GT strength over the low-lying intermediate states (AGTS) 
without the change of their total strength. This latter is mainly determined by the difference of the 
sum rule 3e 2 q (N — Z) and the strength exhausted by the GTS. Nevertheless, because the 2vf3f3 
decay amplitude is formed at the expence of a rather small number of the AGTS (this statement 
follows from the analysis of numerical calculations of Mqt), the consideration for the quasiparticle 
interaction can result in some change of Mqt- The further analysis will allow to specify these reasons 
and, probably, to improve the description of experimental data. 
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Table 1. Explanations are given in the text. 



parent 
nucleus 



type of 
decay 



years -1 m 2 



M GT 
mj 1 



<>GT 



rpexp. 
1 l/2 

years 



rpcalc. 
1 l/2 

years 



76 Ge 
78 Kr 

82 Se 
96 Zr 
m Ru 

100 Mo 
W6 Cd 

U6 Cd 
124 Xe 

i28 Te 
i30 Te 

130 £?a 

136 Xe 
136 Ce 

150 Nd 



P~P~ 



1.317 x 10 



-19 



m 



ec ec 


1.957 


X 


10" 


-21 


(3 + ec 


1.174 


X 


10" 


-21 


p (3 


4.393 


X 


10" 


-18 


pr fr- 


1.953 


X 


10" 


-17 


ee ec 


6.936 


X 


10" 


-21 


P+ec 


1.148 


X 


10" 


-21 


P P 


9.553 


X 


10" 


-18 


ec ec 


1.573 


X 


10" 


-20 


P + ec 


1.970 


X 


10" 


-21 


P P 


8.000 


X 


10" 


-18 


ec ec 


5.101 


X 


10" 


-20 


P+ec 


4.353 


X 


10" 


-21 


p-p- 


8.624 


X 


io- 


-22 


p-p- 


4.849 


X 


10" 


-18 


ec ec 


4.134 


X 


10" 


-20 


P+ec 


1.387 


X 


10" 


-21 


p-p- 


4.870 


X 


10" 


-18 


ec ec 


3.988 


X 


10" 


-20 


p+ec 


6.399 


X 


10" 


-22 


p-p- 


1.200 


X 


io- 


-16 
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m 




















n 







0.0919 
0.0737 



0.0459 
0.0362 



0.0954 



0.0754 



0.0123 
0.0087 



0.0299 



0.0221 
0.0304 



0.0387 

0.0285 

0.0295 
0.0678 
0.1005 

0.1606 
0.1947 

0.0788 
0.0528 

0.0529 
0.0468 
0.0568 

0.0341 
0.0512 

0.0642 



0.131 

0.090 

0.093 
0.324 
0.338 

0.329 
0.319 

0.258 
0.085 

0.090 
0.085 
0.082 

0.088 
0.081 

0.182 



0.9 x 10 21 
1.43 x 10 21 



1.08 x 10 20 

3.9 x 10 19 

> 6.7 x 10 16 
1.15 x 10 19 

> 6.6 x 10 18 
2.25 x 10 19 



7.7 x 10 24 
2.7 x 10 21 

> 4 x 10 21 

> 4 x 10 21 
> 2.3 x 10 20 



1.7 x 10 19 
9 x 10 18 














22] 





[24] 

[24] 







mi 



5.0 x 10 21 

6.2 x 

1.0 x 
2.6 x 

1.1 x 
1.4 x 

8.6 x 

4.1 x 

1.7 x 

1.3 x 

1.2 x 

7.0 x 
8.2 x 

4.1 x 

9.4 x 

7.5 x 

2.2 x 
1.1 x 

9.6 x 
6.0 x 
2.0 x 
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